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O Abstract 

This article reviews the current status of theoretical modeling of electric dipole 
£N| radiation from spinning dust grains. The fundamentally simple problem of dust grain 

rotation appeals to a rich set of concepts of classical and quantum physics, owing the 
i— —i the diversity of processes involved. Rotational excitation and damping rates through 

various mechanisms are discussed, as well as methods of computing the grain angular 
momentum distribution function. Assumptions on grain properties are reviewed. The 
robustness of theoretical predictions now seems mostly limited by the uncertainties 
Q_i regarding the grains themselves, namely their abundance, dipole moments, size and 

q shape distribution. 

iHj 1 Introduction 

Rotational radiation from small grains in the interstellar medium (ISM) has been suggested 
as a source of radio emission several decades ago already. The basic idea was first introduced 
^ by Erickson (1957) [8j and then revisited by Hoyle and Wickramasinghe (1970) [12J and 

Ferrara and Dettmar (1994) Rouan et al. (1992) |23j were the first to provide a thorough 
t— i description of the physics of rotation of polycyclic aromatic hydrocarbons (PAHs), although 

^ not including all gas processes. 

Shortly after the discovery of the anomalous dust-correlated microwave emission (AME) 
j> in the galaxy by Leitch et al (1997) [18], Draine and Lazarian (1998, hereafter DL98) [5j [6] 

suggested that spinning dust radiation might be responsible for the AME, and provided an 
$_i in-depth theoretical description of the process. 

Understanding the spinning dust spectrum in as much detail as possible is important. 
First, the AME constitues a foreground emission to cosmic microwave background (CMB) 
radiation. Second, it provides a window into the properties of small grains, which play crucial 
roles for the physics and chemistry of the ISM. 

Motivated by these considerations and the accumulating observational evidence for dif- 
fuse and localized AME, several groups have since then revisited and refined the DL98 model 
[2J EU QUI EH |29] . New physical processes were accounted for, that can significantly affect 
the predicted spectrum. A publicly available code to evaluate spinning dust emissivities 
(SpDust) is now available, including most (but not all thus far) processes recently investi- 
gate^ 



1 SpDust is available at |http: / / www.sns.ias.edu/ ^yacine / spdust / spdust .html 
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The purpose of this article is to provide an overview of the physics involved in modeling 
spinning dust spectra. We attempt to provide a comprehensive description of the problem 
at the formal level, and let the interested reader learn about the details in the various works 
that deal with the subject. 

This article is organized as follows: Section [2] reviews the basic process of electric dipole 
radiation and the resulting emissivity. We then describe the assumed properties of the small 
grains which are believed to be the source of the spinning dust radiation in Section [3j Section 
[4] discusses the rotational configuration of small grains stochastically heated by ultraviolet 
(UV) photons. Section [5] describes the methods to obtain the distribution of grain angular 
momentum, as well as the various physical processes that affect it. We conclude and mention 
potential future research directions in Section [6] 



2 Basic process 

2.1 Electric dipole radiation of a spinning grain 

Consider a grain with permanent electric dipole moment fi(t) rotating classically, i.e. such 
that its angular momentum is much larger than h (this was shown to be indeed the case even 
for the smallest grains, for which J = L/h ~ 70 |6j). The instantaneous power radiated is 
given by 

p(t) = ^ (i) 

Averaging this power over many rotation cycles, we get 

<^> = ^) = ^E^ 2 ) = ^ £ / = / ^(2™r ;>>>), ( 2 ) 

i i i 

where S^u) is the power spectrum of the i-th cartesian component of fi. From Eq. ^ we 
can directly read off the power radiated per unit frequency: 

— (u\L,^) = —(27rurJ2S^)- (3) 

i 

In general the power spectrum of the electric dipole moment depends not only on the total 
angular momentum of the grain L but also on the the orientation of the grain axes and 
dipole moment with respect to L, which we formally represent by the set of angles vj. 



2.1.1 Spherical grain 

We first consider the simplest case of a freely rotating spherical grain, with isotropic moment 
of inertia tensor 1^ = 18{j. In that case the angular velocity uj$ = L/I is a constant vector, 
which we take along the z-axis, and the power spectra of the components of the dipole 
moment are 

1 

2' 

where /j, ± is the component of /j, perpendicular to u> and z/ = u /(2n) = L/(2tiI) is the 
frequency of rotation. The power radiated in this simple case is then 

dP 2n? 



Sp, = 0, Sf, x = = -H 2 ± 5(v - u ), (4) 



dv 3c 3 



(2vrz/) 4 5(z/ - L/(2tt/)) (5) 
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Figure 1: Euler angles used for the description of an axisymmetric grain. Figure reproduced 
from Ref. 



2.1.2 Axisymmetric grain 

Here we consider an oblate axisymmetric grain with moments of inertia I3 > I2 — I\. We 
describe the orientation of the grain principal axes with respect to its angular momentum L 
with the three Euler angles <fr, 6, i/j pictured in Fig. [TJ 

Between two discrete events that change its angular momentum, the grain can be con- 
sidered as freely rotating. During these periods, the Euler angles change according to 

6 = constant, (6) 

i = y> W 

i, = -(|-£)cos0. (8) 
Electromagnetic radiation is now emitted at the four frequencies </>/(27r), |'0|/(2tt), (<fi ± 

v>)/(2tt) mm- 



dP 



^ 4 sin 2 #5 \v- — \ +^ 4 sin 2 fl5l v- 



dv 3c 3 r I 2tt I 3c 3 r V 2tt 

+ E ^(^^ (1±CO s^L-%«V (9) 



6c 3 1 \ 2ir 

2.1.3 Triaxial grain 

The case of a triaxial grain is unfortunately not analytic. In that case power is radiated 
at a countably infinite number of frequencies. Hoang et al (2011) [H] obtained the power 
spectrum of a freely rotating grain numerically. They show that for a classically rotating 
grain, only a few modes are dominant. 
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2.2 Emissivity 

The quantity of interest to us is the emissivity j u (with units of power per frequency interval 
per unit volume per steradian), which we obtain from integrating the power over grain size 
and shape distribution dn gv /da (where a is meant to formally represent the characteristic 
size and three-dimensional shape of a grain), as well as the electric dipole distribution^] 
/(/i|a), convolved with the probability distribution for the angular momentum and rotational 
configuration f(L, zu\a; /j,). In the absence of any preferred direction, the dependence on L 
is in fact only through its magnitude L. This is the case if there are no magnetic fields and 
anisotropic radiation fields, or if none of them are efficient at aligning the grains. In what 
follows we will assume perfect isotropy of space. The emissivity is then given by 

1 f dn f f dP 

3u = 4n / da ~da / ^^f^a) J d z L dzu f(L,zu\a, ^) — (u\L,zu,a, (10) 

Finally, we note that the actual observable is the radio intensity, given by the emissivity 
integrated along the line of sight, 

I v {h) = J j u (sh)ds. (11) 

At each point in space, the emissivity depends on the local environmental conditions (density, 
temperature, ionization degree, ambient radiation field, as well as grain abundance), as we 
shall see below. Therefore predicting the spinning dust spectrum along a line of sight requires 
modeling the ISM properties (see for example Ref. [28J ) . Ref. pj] evaluate the effect of 
turbulence on the spinning dust spectrum, and find that the effective emissivity (averaged 
over the probability distribution for the compression factor along the line of sight) can be 
shifted to larger frequencies and enhanced by several tens of percent. 

We do not deal with this aspect in this review as it does not belong, per se, to the field 
of spinning dust theory, but rather to the larger field of ISM modeling. It is however crucial 
to accurately model the environmental spatial variations in order to get precise predictions. 



3 Grain properties 

3.1 Abundance and size distribution 

The small grain abundance is determined primarily from observations of the wavelength- 
dependent extinction^] and the 3 — 25 /mi emission, attributed to various vibrational modes 
of PAHs (for a review on interstellar PAHs and their properties, see for example Tielens 
(2008) [26]). Observations require a few percent of the interstellar carbon to be locked in 
PAHs, and a significant population of very small grains (less than ~ 1 nm in size) to reproduce 
the strength of the 3-12 urn. bands, which are emitted by small grains stochastically heated 
to large temperatures [7]. Following DL98, Li and Draine (2001) [H] and Weingartner 

2 We use repeatedly the letter / to denote probability distribution functions. The arguments of / should 
always make its meaning unambiguous: f(X, Y\a,(3, . . .) denotes the differential probability distribution for 
the variables X, Y given the values a, /3, . . . of some other variables. 

3 The UV extinction indicates the presence of nanodust but does not give any detailed information on 
grain sizes. Only the IR emission allows to constrain the small grain size distribution, as discussed in Li & 
Mann (2012) [3T]. 
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and Draine (2001) [27] have adopted a log-normal distribution in grain radius, starting and 
centered at a min = a = 3.5 A (corresponding roughly to a coronene molecule with formula 
C24H 12 ) and variance a = 0.4 in log a, and showed that such a distribution reproduced the 
infrared emission well. Observations tend to indicate that PAHS are less abundant in dense 
clouds than in the diffuse ISM. 

Note that we assume that the smallest grains in the ISM are mostly PAHs, but a popu- 
lation of ultrasmall silicate grains is not completely ruled out by observations [20] . 

3.2 Shapes 

PAHs may take a variety of shapes, from disk-like to nearly linear. They are not necessarily 
planar: for example if one of the hexagonal carbon rings is replaced by a pentagonal ring, 
they are bent and become three-dimensional. Above a certain size, PAHs may form irregular 
clusters, and eventually, large three-dimensional grains. 

The exact distribution of shapes is largely unknown. The lowest-frequency IR emission 
bands in principle carry information about the individual grains, and seem to indicate that 
PAHs may be dominated by a few well-defined molecular structures, although not conclu- 
sively [2"o] . 

The smallest grains dominate the spinning dust spectrum (they can be spun up to larger 
frequencies and hence emit more power). It is commonly assumed that these grains are 
nearly planar up to a spherical-equivalent radius a = 6 A, corresponding to 100 carbon 
atoms. The peak of the spinning dust spectrum is not very sensitive to the exact cutoff 
between planar and spherical grains. 

3.3 Permanent dipole moments 

In principle a consistent prescription should be given for small grains, that gives the precise 
nature of the grain, hence its shape (or rotational constants) and permanent electric dipole 
moment, which can be computed quantum-mechanically for small enough molecules. Such 
computations were carried out by Hudgins, Bauschlicher and Allamandola (2005) [13] for 
nitrogen-substituted PAHs. They found typical permanent dipole moments of a few Debyes, 
depending on the precise position of the substituted nitrogen atom. 

Eventually, observations will hopefully allow for a more precise determination of the pop- 
ulation of PAHs and their properties. We are currently far from having a definite handle on 
such refined properties of small grains, and an empirical distribution of dipole moments is 
required. Following DL98, more recent models assume a three-dimensional gaussian distri- 
bution of dipole moments, with variance 

(ti 2 ) = N &t (3 2 + (eaeZ) 2 , (12) 

where /3 ~ 0.4 Debye and e ~ 0.01. The first term, largely dominant, accounts for the 
permanent dipole moment and the second term accounts for charge displacement in ionized 
grains. We repeat that this distribution is largely ad hoc and may be far from reflecting 
reality, except (hopefully) for the characteristic permanent dipole moment. 
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4 Rotational configuration 

4.1 Fast vibration-rotation energy transfer 

In principle one should solve for the distribution of angular momentum and rotational con- 
figuration f(L, vj) at once. Indeed, a priori, the same processes that change the angular 
momentum may also change the relative orientation of the grain at a similar rate. Grains can 
exchange angular momentum with several "baths" , all characterized by different characteris- 
tic "temperatures" , and the resulting overall distribution function cannot easily be factored 
into a pure angular momentum part and a pure rotational configuration part (one can al- 
ways formally write the factorization, f(L,vj) = f(L)f(vj\L), but one cannot in principle 
compute the two factors independently). 

The situation is much simplified if one single process is very efficient at changing the 
rotational configuration, on timescales much shorter than the overall timescale to change the 
angular momentum. If this process is characterized by an equilibrium temperature T OT , then 
one can indeed compute the probability distribution f(vj\L) independently: 



f(zu\L) oc exp 



E ro t(L, vj) 



kT, 



VJ 



(13) 



where the rotational energy is most easily written using the projections of the angular mo- 
mentum on the grain's principal axes with principal moments of inertia Jj, 

i i 

where we have defined U = Li/L, the normalized projection of the angular momentum 
along the axis i. One can formally identify vj <h- {/«} (at least in some time-averaged 



sense), and Eq. (13) then uniquely determines the probability distribution for the rotational 
configuration, given a value for the total angular momentum L. The angular momentum 
distribution f(L) can then be obtained from transition rates averaged over the rotational 
configuration with the known distribution f(vj\L). 

Luckily Nature does provide us with such an efficient process to change vj at constant 
L: internal vibrational-rotational energy transfer (IVRET) (see for example [221 EH 122] an d 
[2HHI] f° r application to spinning dust modeling). The first detailed studies of how internal 
relaxation may affect grain alignment were carried in Refs. [T5l \T7\ . 

Following the absorption of an ultraviolet (UV) photon, small grains get heated up to 
large vibrational temperatures T v ;b (the notion of temperature is not well defined for the 
smallest grains, so we mean temperature as a characteristic energy per degree of freedom). 
IVRET leads to a rapid energy exchange between vibrational and rotational degrees of 
freedom, at constant angular momentum, so that during a thermal spike, the distribution 



f(vj\L) is given by Eq. (13) with T OT = T v ;b. As the grain cools down by emitting infrared 
photons, its vibrational temperature decreases, until the grain reaches its fundamental vi- 
brational mode, with typical energy E /k ~ 100 K for the smallest grains. IVRET is only 
active as long as the density of vibrational states is large enough that there exist transitions 
at frequencies near the rotation frequency. Therefore energy exchange freezes at a character- 
istic temperature T ro = Tf r > E /k, probably of at least a couple hundred Kelvins, and this 
temperature characterizes the final rotational configuration following a thermal spike. For 
this distribution to remain valid at all times, it is necessary that the rate of absorption of UV 
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photons is larger than the rate of change of angular momentum. This is indeed the case in 
diffuse environments (see Table 1 of Ref. [21]), where absorption of UV photons is a few times 
faster than angular momentum changes; this difference of timescales gets more pronounced 
as the ambient radiation field increases. In dense and under-illuminated clouds, however, 
the rate of absorption of photons is not large enough to maintain the distribution (13), and 
the rotational configuration needs in principle to be computed from the full f(L, vj). 

Let us now discuss the implications of Eq. (13). The angular momentum itself is has a 
characteristic value L pea k, and therefore the rotational energy is of order E rot ~ L^ cak /I = 
kT^. If T OT <C Tl, the most probable rotational configuration will be the one minimizing 
the energy, i.e. where the grain rotates about its axis of greatest inertia. In the opposite 
case where T ro ^> Tl, all rotational configurations become equiprobable (of course when 
converting to actual angles one needs to be careful of using the appropriate phase-space 
volume dvj). 

For example, an axisymmetric oblate grain with J 3 > J 2 = I\ has a rotational energy 



L 



2 



cos 2 6 



(15) 



h 

h, 

where 9 is the angle between L and the axis of greatest inertia. If <C Tl, the most 
probable configurations are 9 = or 7r. In the case where T ro 3> Tl, we obtain that 
f(9\L) oc sin# (here we used dvj = d<fi sin 9 d9 in the usual spherical polar coordinates). 

SpDust only allows for the two limiting regimes T^ — > and T ro — > oo, bracketing the 
range of possibilities. The authors of Ref. [TT] explore the effect of continuously varying 
T OT , interpolating continuously between the two extreme regimes. In general there will be 
a different temperature T^ for each grain size and depending on the environment, but the 
precise modeling of this parameter has not been addressed in the literature yet. In what 
follows we shall only discuss the two limiting regimes. 



4.2 Implication for the emitted power at fixed angular momentum 



The last integral in Eq. (10) can be rewritten as 

dP 

d z Ldvjf[L,vj) — [L,-uj) - 



(16) 



(17) 



where the power averaged over the rotational configuration is 

/dP> 
\~du 

In the case of a grain rotating about its axis of greatest inertia /3, the averaged power 
collapses to 

( d ^) = ^@xv)Vj{v - L/(2nh)), T OT -> 0, (18) 

which is identical to the case of a spherical grain. If we now consider an oblate axisymmetric 
grain with T ro — > oo, we obtain, using the results of Section 2.1.2, and averaging over 
isotropically distributed nutation angles 9, 



f dP 
(L) = j dvjf{vj\L) — {L,vj). 



(-) 
\du/ 



^(2vr, 1 ) 4 ^-^) + §(2vr,r 



1 - 



v 



^13 



"13 



~\-v>,<v<v\+v\?, 
2^13 



T-rr, -)■ OO, 



(19) 
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where the function 1 is unity where its subscript is valid and zero elsewhere, and we have 
defined the two frequencies 

« s 2J7? (20) 

In the case of a planar grain with with J 3 = 21 \ and v\ = 2u 3 = 2z/ 13 , the power radiated by 
a wobbling grain is emitted at characteristic frequencies about twice as large as in the case 
of a grain rotating primarily about its axis of greatest inertia. The integrated power is, in 
the former case (and assuming 1/3 = ^13) 

P(T ro 00) = ^^i) 4 + 10^(2™ 3 ) 4 , (22) 

which is about 10 times larger, at equal angular momentum, than the power radiated by a 
grain rotating mostly about its axis of greatest inertia, if jjr\ = 2//n. 

One must not forget, however, that the angular momentum distribution itself depends 
upon the rotational configuration f(zu\L), since one must use this distribution to average 
transition rates. We shall see in the next section that the effect of randomized rotational 
configuration is to lower the characteristic angular momentum L. 



5 Angular momentum distribution 

To determine the angular momentum distribution, we need to evaluate the differential transi- 
tion rates T(L — > V) between different values of the angular momentum magnitude, defined 
such that T(L — > L')AL' is the rate of transition from an initial angular momentum L to 
a final angular momentum in the interval [V , V + AL']. These rates are averaged over the 
rotational configuration for zu discussed above. The steady-state distribution function f(L) 
should then in principle be obtained from the integral master equation, 



df(L) 
dt 



J f(L')T{L' L) — f{L)T(L — >■ L') dll = 0, for all L, (23) 



where we have defined f(L) = 4irL 2 f(L) so that f(L) is the distribution function for the 
magnitude of L (whereas f(L) = f(L) is the distribution function for the vector angular 
momentum, even if it only depends on its magnitude due to isotropy). This equation is, 
clearly, rather cumbersome to solve and below we present a simpler (if approximate) method 
of solution, based on the Fokker-Planck equation. Section [KT] provides a formal introduction 



to the problem, and actual physical mechanisms are discussed in Section 5.2 



5.1 The Fokker-Planck equation 
5.1.1 Derivation 

In general transition rates are not significant for arbitrarily large values of AL = L' — L: there 
always exists some characteristic ALo such that T(L — > L + AL) decreases rapidly for |AL| > 
AL . An important simplification can be made if the scale AL is much smaller than the 
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characteristic scale over which both the distribution function varies and the rates themselves 
vary - the distribution function being unknown a priori, the validity this assumption has in 
principle to be checked a posteriori. If this is the case, we may expand the first term in the 



integral of Eq. (23), setting I! = L + AL: 



d 

f(L + AL)T(L + AL — >■ L) « f(L)T(L -»■ L - AL) + AL— f(L)T(L -»■ L - AL) 



+ 



1 d 2 



3L V 

/(L)r(L -> L - AL) 



(24) 



Plugging this expansion back into Eq. (23), we see that the term linear in f(L) cancels out 
with the second term of the integral (the integrand being an odd function of AL). Recalling 
that f(L) = AnL 2 f(L), we finally obtain 



AttL 



2 df(L) 



Of 



d_ 
dL 



d(AL) 
dt 



4nL 2 f(L) 



+ 



1 d 



2 r 



2<9L 2 



d(( AL) 2 ) 
dt 



4nL 2 f(L) 



0, (25) 



where we have defined the rate of angular momentum drift (the opposite of which is the rate 
of angular momentum dissipation or damping), 



d(AL) 
dt 



AL T(L — )• L + AL)d(AL), 



and the rate of angular momentum diffusion (also termed fluctuation or excitation), 

d((AL) 2 } 



dt 



J (AL) 2 r(L -> L + AL)d(AL). 



(26) 



(27) 



Equation (25) is known as the Fokker-Planck equation and has a broad range of applications 



in physics (see for example Chapter 6 of the book [3] by Blandford and Thorne). In the 
context of grain rotation, this approach was used in Refs. [HJ [HI El El], and is the basic 
equation solved in SpDust. Solving this equation is much simpler than solving the full 
master equation, and it may even have a simple analytic solution if the excitation and 
damping rates are simple enough. 



5.1.2 General form of the rates and solution 

General processes besides electric dipole radiation damping. Most processes through 
which grains may change angular momentum (except for electric dipole radiation itself, to 
which we shall come back later on) are characterized by a damping timescale r such that 



d(AL) 
dt 



L 

T 



and have an isotropic and constant diffusion rate of the form 

d(ALiALj) 



dt 



Taylor-expanding AL to second order in AL, we obtain 

(AL) 2 — (AL ■ L) 2 



AL= IL + ALI - L« AL-L 



2L 



(2f 



(29) 



(30) 
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The excitation and damping rates for the magnitude of the angular momentum are therefore 



d(AL) 
dt 

dt 



oi/3 - L 2 



Lt 



3r' 



(31) 
(32) 



The form of these coefficients stems from the fact that only longitudinal excitation generates 
a true diffusion in the magnitude of L, whereas excitations perpendicular to L lead to a 
systematic increase of the magnitude of the angular momentum, and hence appear as a 



positive drift rate in Eq. (31) 



If only one single process was interacting with the grains, their steady-state distribution 
would then be the Maxwellian with three-dimensional variance (L 2 ) 



f{L) oc exp 



3L 2 



(J, 



(33) 



as can be seen from inserting the rates (31) and (32) into the Fokker-Planck equation (25). 



Often, however, there is not a single process that dominates both excitation and damping. 
Since transition rates add up linearly for independent processes, so do excitation and damping 
rates. If all rates are of the form (28), (29), then the final distribution is still Maxwellian, 



with a variance weighted by the characteristic rates of the various processes indexed by a: 



f(L) oc exp 



3L 



(Jl 



21 



y 



-1^2 



_-l 



.-1 



(34) 



Damping through electric dipole radiation. One process behaves differently from 



Eqs. (28), (29): the damping of angular momentum through electric dipole radiation itself. 
Since the rotational energy is proportional to L 2 and the radiated power is proportional to 
L 4 , the rate of angular momentum damping scales as L 3 . We may write it in the form 



d(AL) ] 
dt 



led 



L 2 L 

2 : 



(35) 



where we may take the variance a L to be that given by Eq. (34). This defines r oc j as the 



characteristic timescale to damp an angular momentum of order o~l through electric dipole 
radiation (note that this definition of r e d is different from the ones adopted in Refs. [HI [2] , 



where a L = 3IkT gas was specifically used in Eq. (35) to define r e d). Every damping process 



has in general an associated excitation process, and vice versa. In the case of electric dipole 
radiation, the associated fluctuation in angular momentum is due to absorption of and decays 
stimulated by microwave photons (dominated by Cosmic Microwave Background (CMB) 
photons in the diffuse ISM). To our knowledge this process was only considered in Ref. [1] 
and we shall get back to it in the next section. 

Accounting for the damping only for now and including this additional damping into the 
Fokker-Planck equation, we obtain the solution 



f(L) oc exp 



L 2 



3 r L 4 

4 



a 2 2 r cd erf 



(36) 
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The most likely angular momentum for this distribution is such that the net drift rate 
vanishes, 

-kpeak = -7T / (37) 

Oil h i 4 t 



1 + V 1 + 3^ 

For r e d *C r, which is often the case for the smallest grains, the most likely angular mo- 
mentum results from equilibrium between damping through electric dipole radiation and 
excitations through other mechanisms, and is approximately 

Lle, k ^^°l ^ed«r. (38) 

In that case, the characteristic damping time at the peak is r rot = y/rr ed . In general, one 
can define a characteristic damping time as 

r rot min[r, y/rr ed \. (39) 

We can rewrite the peak angular momentum in terms of r rot as 

LU ~ T -f°l (40) 

5.1.3 Limitation of the Fokker-Planck approach: impulsive torques 

The Fokker-Planck equation is a diffusion equation, and its validity is limited to processes 
that change the angular momentum by small increments. In this section we formally discuss 
in which cases it may break down. 

Let us now consider some stochastic interaction process a (in practice, collisions with 
passing ions [TO] ) with rate r~ u , such that the characteristic angular momentum exchanged 
with a grain at each interaction has variance 

(ALAL^a = (41) 

T~n O 



The diffusion rate for such an interaction is indeed, formally equal to that of Eq. (29). 
However, this process can only be considered as diffusive if ((AL) 2 ) a <C L 2 eak . This translates 
to the condition 

r -i tT 2 

7-coii < r rot _ 1 2 L ■ (42) 

T a a L,a 

If the process a is the dominant excitation mechanism, t~ x o\ = t~ 1 <j 2 l a and the condition 
for the validity of the diffusion approximation is that r co ii T rot , the characteristic time to 
change the angular momentum. 

The issue of impulsive torques was addressed by Hoang et al [TU]. Instead of solving a 
Fokker-Planck equation, Hoang et al start with the Langevin equation, of the form 



dL= ^ dt+ .mm dq , ,43) 

dt V dt 



where dq is a random variable with variance ((dq) 2 ) = dt. They then solve it numerically 
to obtain L(tj) at a set of discret timesteps t{. The distribution function f(L) is then 
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obtained from the histogram of values of Lj after a long enough evolution. In the form (43), 
the Langevin equation is exactly equivalent to the Fokker-Planck equation, in the sense 
that it assumes infinitesimal torques. However, it is simple to generalize this treatment 
to include impulsive toques. First, one draws the interval between two collisions from the 
Poisson distribution with mean t co \\. Second, an angular momentum AL, is drawn from 
an isotropic distribution with variance ((AL) 2 oll ). This method allows to include random 
impulsive torques in addition to the quasi-continuous torques. We defer the discussion of 



Hoang et al's results to Section 5.4 



5.2 Excitation and damping rates for various mechanisms 

In this section we describe the principal mechanisms that excite and damp the grains' rota- 
tion. Since the detailed calculations are already worked out in various papers [21 1231 El I2"4] . 
here we limit ourselves to giving a semi- qualitative description of each process and order of 
magnitude estimates for the relevant rates. Since the smallest grains are producing the peak 
of the spinning dust spectrum, all numerical evaluations are normalized to the characteristic 
radius of coronene, a ~ 3.5 A. 



5.2.1 Collisions 

Collisional interactions of grains with gas atoms, molecules or ions is perhaps the most 
intuitive of angular momentum transfer processes, even though the microphysical details 
could be very complex (see for example the discussion in Section 4.2 of Ref. [23J ) . Impactors 
with density n imp reach the grain with a rate 

T coll ~ ^impTT&max^inj (44) 

where vr6^ ax is the effective collisional cross-section (6 max being the maximal impact param- 
eter for which a collision occurs) and u m is the characteristic velocity of the impactors at 
infinity. As they impact the grain, they stick to its surface, providing there are available 
adsorption sites. The random angular momentum transferred to the grain has variance 

((AL) 2 ) 

in ~ ™imp"max Vin) , (45) 

where m; mp is the mass of the impactor. 

The attached impactors are ejected from the grain's surface following the absorption of 
UV photons that heat up small grains to large temperatures (this is the process of photo- 
evaporation). Because ions are in general more electronegative than large molecules, they 
leave the grain surface as neutral species. Here again, they give a random recoil to the grain, 
leading to 

((AL) 2 ) out ~ (mimpGWout) 2 , (46) 

where v out is the characteristic velocity of ejection, related to the dust grain temperature 
T vib following thermal spikes by m imp v^ ut ~ 

In addition to a random component, ejected particles systematically decrease the angular 
momentum of the grain: if their ejection velocity is random in the rotating grain's frame, 
they carry on average an angular momentum 

2 

at 2 m impO' T m \mp T f An\ 

AL out = m imp a u ~ — L ~ L, (47) 

1 ^grain 
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where u is the rotation rate of the grain and / ~ rn gra i n a 2 is its characteristic moment of 
inertia. In steady-state, the rate of ejections equals the rate of collisions and therefore 

<*< Ai >- -r-l^L. (48) 



Lib //fc grain 



From this expression we see that the characteristic timescale for damping the angular mo- 
mentum through ejection of colliding gas particles is 



^grain //in\ 
T ~ — — T coU . (49) 



Using the definition (32), we see that the characteristic variance in angular momentum that 



would stem from collisions alone is 

r / 

<y\ ~ ((A£) 2 ) ~ 2 [( m imp&maxWin) 2 + (™imp<Wout) 2 ] 

7"coll ^imp^ 



2 

max 



(I 



2 



kT\ n -\- kT c 



out 



• (50) 



We see that collisions tend to drive the angular momentum distribution to a thermal distri- 
bution with temperature 

1 



- L,coll 2 



'-,2 



max rp , rp 
~ -L in i -L o 



(I 



2 



(51) 



The maximum impact parameter (or effective cross-section) for incoming particles depends 
upon the charge state of the grain and impactor. It can easily be determined at an order-of 
magnitude for a spherically symmetric interaction potential V(r). Requiring energy and 
angular momentum conservation, we find 



^ 2 

'max 



..m^.m, (52) 

In the case of a repulsive interaction V > 0, this should be understood as 6 max ~ a if 
V(a) kT in and in the opposite case. 

• For neutral impactors and neutral grains, 6 max = a. 

• For neutral impactors and charged grains (with charge Z g e), the collisional cross-section 
can be determined from equating the kinetic energy of the incoming particle to the potential 
energy of the attractive induced-dipole interaction, 

1 Z„ 2 e 2 

V(r) = -2«-£-> (53) 

where a is the polarizability of the impactor. Typically, a ~ 1 A 3 so the focusing factor 
(frmax/ a ) 2 is of order ~ 1 + (T/10 3 K) -1 for a w 3.5 A. This potential must also be accounted 
for when evaluating the escape probability of neutral particles ejected from ionized grains. 

• For positively charged impacting ions, the dominant interaction is the attractive Coulomb 
attraction with negatively charged grain^J 

I 7 7 \p 2 

V( r ) = - 1 1 gl , (54) 
r 



4 Whenever collisions with ions are relevant, a significant fraction of grains are negatively charged by 
colliding electrons, so the cation-PAH - collisions are in general dominant over cation-PAH° and cation- 
PAH+. 
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corresponding to a large focusing factor 

2 I7.7l„2 



frmax \ \ZiZ„\e _„ / a \ I T 



a 



akT V 3.5 A/ V10 3 K 



50 — x 77^7 • (55) 



Therefore we see that collisions with ions may overcome collisions with neutrals even for 
relatively small ionization degrees. 

5.2.2 Plasma excitation and drag 

Ions can exchange angular momentum with the grains at a distance, without necessarily 
colliding with them, by exerting a torque on their permanent electric dipole moment, 

dL _ Zie 

The characteristic interaction timescale being At ~ b/v, the variance of the angular momen- 
tum change for each interaction event is 

« Ai)2> ~(^) 2 - (57) 

Integrating over impact parameters b and velocities of ions v, one obtains the rate of angular 

momentum diffusion 

<(AL) 2 > 



/Ltx.E~//-^-. (56) 



(H n^eyj^lnA, (58) 

where the order-unity Coulomb logarithm In A = ln(6hi/&io) appears because of the logarith- 
mic divergence of the integral over impact parameters. In practice the lower cutoff must be 
set to the maximal impact parameter leading to a collision (denoted 6 max in the previous 
section), in order not to double-count the angular momentum transfer in collisions. The 
upper cutoff comes from the fact that when the interaction timescale b/v is much larger 
than the rotation timescale u , the torque along the angular momentum vector averages 
out to zero j2]. 

To obtain the damping rate from first principles would require accurate evaluations of the 
back-reaction of the grain on the ions' trajectories, leading to a small asymmetry between 
trajectories increasing the magnitude of L and those decreasing it. However, a powerful 
theorem, the fluctuation- dissipation theorem [31 Sj, allows us to very simply evaluate the 
dissipation rate from the fluctuation rate if the interaction is with a thermal bath. Put 
simply, excitation and damping must balance in such a way that, if only the thermal process 
considered was at play, the resulting distribution would also be thermal, of the form f(L) oc 
exp[— E(L)/kT] with the same temperature T as the bath. 

In the case of a grain rotating about its axis of greatest inertia, and using the notation 



of Section 5.1.2, the damping timescale for plasma drag r p must therefore be such that 



3r p ^^ = al p = QIkT gas . (59) 

The case of a wobbling grain is a little more complex, since the temperature for the rotational 
configuration is set by the internal relaxation process and is in general different from the gas 
temperature. However, one can use a closely related principle, that of detailed balance, to 
compute the proper damping rate given the tensorial excitation rate. Details can be found 
in Ref. 
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5.2.3 Emission of infrared photons 

Every time a small dust grain absorbs a UV photon, it gets into a highly excited vibrational 
state from which it decays by emitting a cascade of infrared (IR) photons, typically about a 
hundred per absorbed UV photon. Each one of the emitted IR photons carries one quantum 
of angular momentum, so its angular momentum squared is L 2 = 2h 2 . If photons are emitted 
isotropically with a energy flux F u , the rate of diffusion of angular momentum is then 

d((AL) 2 ) 1 - f , F v 
u ' 1 - -2h 2 / (60) 



dt 3 J hv 

A ro-vibrational transition from the excited state (v, J) [v denoting the vibrational configu- 
ration and J the rotational quantum number] to an other state (V, J + A J) has a transition 
frequency 

u = u - ^AJ, (61) 

where v§ is the transition frequency for (v, J) — > (v', J) and u = hJ/I is the angular rotation 
frequency. Quantum-mechanical transition rates are proportional to the transition frequency 
cubed. Summing over the three allowed transitions A J = 0, ±1 (and assuming they have 
nearly the same matrix elements), the net rate of angular momentum drift relates to the 
rate of photon emission iV 7 through 

J = (62 > 

which implies, in the case of a grain rotating about its axis of greatest inertia, 

dt (2n) 2 I J v 2 y ' 

A classical calculation can be found in Refs. |6] (with missing factors of two for both damping 
and excitation rates), [2] (with a missing factor of two for the excitation rate) and [21]. A 
fully rigorous quantum-mechanical treatment can be found in Refs. [21 H [231 [2H]- They are 
perfectly equivalent since ISM grains are classical rotators. 



The main difficulty in correctly evaluating Eqs. (60) and (63) is that one must be able to 
compute the infrared spectrum with high accuracy, especially at long wavelengths (where it 
is not well constrained by observations). 

5.2.4 Electric dipole radiation and absorption of CMB photons 

Damping rate A grain emitting electric dipole radiation also radiates away angular mo- 
mentum. Classically, the radiation reaction torque is given by 

d(AL). 2 .. ... 

-V 1 - = -s> x (64) 

where the averaging is over the quasi-periodic rotation of the grain. For a grain rotating 
about its axis of greatest inertia, the result is 



d(AL) 2/4 / L x ' 



dt Ied 3c 3 V/; 



'3 



L, T ro 0. (65) 
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In the case of a wobbling disk-like grain with completely randomized nutation state, the 
corresponding damping rate is (assuming I 3 = H\ for simplicity) 



d(AL) 
dt 



led 



82/xj 32/jjj 
45c 3 + 9c 3 



T„ -)■ oc. 



(66) 



The characteristic timescale to damp an angular momentum L = cr^ [defined in Eq. (35)] is 
therefore such that 



'cd 



3c 3 Jf ' 
82/4 + 32/ijj 



r 



o. 



45c 3 9c 3 



°1. 



T_ -)■ oo. 



(67) 
(68) 



We see that the electric dipole radiation damping timescale is typically about 5 times shorter 
if the grain is wobbling than when it is rotating primarily about its axis of greatest inertia. 



Excitation rate The associated excitation mechanism comes from the absorption of CMB 
photons pp. For simplicity, here we only consider the case of a grain rotating abouts its axis 
of greatest inertia. 

Quantum mechanically, the damping of the angular momentum is due to spontaneous 
decays J — > J — 1, with rate Aj t j_i, so we may rewrite 



d{AL) 
dt 



led 



-HAjrxL. 



(69) 



In addition to these spontaneous decays, stimulated decays take place, as well as absorptions 
of CMB photons. For large J, the net angular momentum change due to these transition 
nearly cancel out (so the drift is essentially due to spontaneous decays). However, the rate 
of excitation parallel to the angular momentum axis is (in the limit that J> 1), 



d((AL l{ ) 2 } 



dt 



led 



2h 2 n~ / (i>)Ajj_i = 2hn 1 {i>) 



d(AL) 



dt 



led 



where v = L/{2ttI 3 ) is the transition frequency and 

, s 1 



(70) 



(71) 



is the photon occupation number at the transition frequency. The CMB temperature is 
T 7 2.73 K, corresponding to a frequency kT^/h pa 57 GHz, so the photon occupation 
number is of order unity at characteristic grain rotation frequencies of a few tens of GHz. 
For a characteristic angular momentum L pea k ~ <Jl(t cc \/t) 1 I a (see discussion in Section 



5.1.2), the ratio of excitations by CMB photons to other excitations is of order 

d((AL)2) 



dt 



led 



d((AL)2) I 
dt l tot 



L 



L 



peak 



n 7 (i/). 



(72) 



For a coronene grain rotating at 30 GHz, the rotational quantum number is typically J = 
L/h pa 70. We therefore conclude that excitations by absorptions of and decays stimulated 
by CMB photons are subdominant, having an effect of the order of a few percent, with a 
greater importance in regions where grains are slowly rotating. 
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Table 1: Dominant excitation and damping mechanisms for the smallest grains considered 
(a = 3.5 A), as a function of idealized environment (see definitions in Table 1 of Ref. [B]). Two 
or more mechanisms are written down if they are of comparable importance, by decreasing 
order of importance, "e.d." stands for electric dipole radiation torque. 



Phase 


DC 


MC 


RN 


PDR 


excitation 


coll. (neutrals, ions) 


coll. (ions) 


IR 


coll. (neutrals) 


damping 


e.d., coll. (neutrals) 


plasma drag 


e.d., IR 


e.d., IR, coll. (neutrals) 



Phase 


CNM 


WNM 


WIM 


excitation 


coll. (ions, neutrals) 


coll. (ions, neutrals), IR 


coll. (ions) 


damping 


e.d. 


e.d 


e.d 



5.2.5 H 2 formation and photoelectric ejection 

Draine and Lazarian j6] considered the random torques exerted on grains as molecular hy- 
drogen is formed on their surface and subsequently ejected, and found that this effect was 
sub dominant. 

Similarly, the rotational excitation due to photoejection of electrons following UV photon 
absorption is a subdominant excitation mechanism. 



5.3 Dominant excitation and damping mechanisms as a function 
of environment 

The relative importance of the various mechanisms described above depends upon the precise 
environmental conditions, i.e. the gas density, temperature, ionization state, and ambient 
radiation field. Note that these parameters also affect the rotational transition rates through 
their dependence on grain charge. Since the timescale for grains to change charge is in general 
shorter than the timescale to change the grain angular momentum (though they are in fact 
comparable for the smallest grains, see Fig. 3 of Ref. [U]), excitation and damping rates must 
be averaged over the grain charge distribution function^} 

We list in Table [T] the dominant excitation and damping mechanisms for the smallest 
grains in the various idealized environments defined in Table 1 of Draine & Lazarian [6] . It 
can be seen that every mechanism discussed above can be dominant under some conditions, 
and several may be of comparable importance in some regions. In diffuse ISM phases, electric 
dipole radiation torque is systematically the dominant damping mechanism, and collisions 
(in general with ions) are almost always the dominant excitation mechanism. 



5.4 Effect of impulsive torques 



We discussed in Section 5.1.3 how to characterize the importance of impulsive torques. In 



this section we discuss specifically the case of the warm ionized medium (WIM), where 



5 As a consequence the electric dipole radiation should not be correlated with indicators of grain charge, 
such as IR line strength ratios. However, since charging time and rotational decay time are comparable 
for the smallest grains, in practice there could be some level of correlation. Quantifying this would require 
solving for f(L, Z) simultaneously, a problem not adressed in the literature. 
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collisions with ions are frequent and the rotational damping time is short. 

The WIM is characterized by a large gas temperature T w 8000 K, an a fully ionized 
gas at low density, nn+ ~ 0.1 cm~ 3 . Collisions with ions provide the dominant excitation 
mechanism. Grains are mostly negatively charged due to the high rate of sticking collisions 
with high-velocity electrons. For a coronene molecule, the characteristic time between ion 
collisions and the characteristic rotational damping time at the peak angular momentum 
Trot = y/TT c d turn out to be comparable^ of order a few years. This indicates that the dif- 
fusion approximation is not strictly correct, and impulsive torques may affect the rotational 
distribution function 

Hoang et al pi)] provided a detailed calculation for the effect of impulsive torques by 
solving a generalized Langevin equation. We reproduce the angular momentum distribution 
function they obtain in Fig. [2} We see that impulsive torques significantly enhance the high- 
frequency tail of the distribution function, due to grains rotating near the peak frequency 
being impulsively spun up to larger rotation rates. This enhancement is mostly unobservable 
because the vibrational emission from large grains dominate at these frequencies. More 
importantly, Hoang et al. found that the peak emissivity is enhanced by about 23% for the 
WIM [and only 11 % for the warm neutral medium (WNM)], although the peak frequency 
remains unchanged. This effect is therefore marginally important for the WIM and should 
be included in precise modeling toola^} 



6 Hoang et al [TO] compare t co ii to t c( j (which is the rotational damping time at L — y/3IkT 3> L pea k and 
in this specific case turns out to be equal to our t 0( j within a factor of a few), and find r cc i <C t co u, by more 
than two orders of magnitude. The correct comparison should be with r rot (the actual rotational damping 
near L pea k), which is in fact comparable to r co ii, explaining the marginal importance of the effect near the 
peak of the distribution function. 

7 This effect is not, as yet, included into SpDust. 
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5.5 Effect of grain wobbling 

A more important effect on the spectrum is that of increasing the characteristic internal 
temperature T OT , which makes the grains wobble rather than simply spin about their axis 
of greatest inertia. It is instructive to make a basic estimate of the effect from simple 
considerations. 

The rotational energy of an axisymmetric grain is given by Eq. (15). Depending on the 
value of T OT , the relation between mean rotation energy (averaged over the distribution of 
nutation angles) and total angular momentum is (assuming I 3 = 2Ii) 

L 2 = 2I 3 (E rot ), T ro ^0, (73) 

L 2 = I x 2I 3 (E TOt ), -> oo, (74) 
5 

If interacting with a bath of characteristic temperature 2~b a th, grains tend to have a charac- 
teristic rotational energy 

(E Iot ) « kT hath , (75) 

regardless of their internal temperature T ro . Therefore, the characteristic angular momentum 
variance (defined in Eq. 32) is 

a|(T ro ^oc)^cx£(T ro ^O). (76) 

The excitation rate being roughly independent of the actual angular momentum, we deduce 
that the damping timescale must scale in a similar fashion as er|, i.e. 

r(T w ^oo)«|r(T w -►()). (77) 
5 

We saw previously that the rate of electric dipole damping is about 5 times larger in the 
wobbling case, at equal angular momentum. The characteristic electric dipole damping 
timescale defined in Eq. (35) is therefore 

r cd (T ro ->■ oo) « ^r cd (T ro ->■ 0). (78) 



Finally, the most likely angular momentum, in the case r e( j <C r, was given in Eq. (38) and 
is therefore such that 

£peak(7~ -> oo) « ^ cak (T OT -> 0). (79) 

The peak frequency is linear in the peak angular momentum, and at equal angular momentum 
it is ~ twice as large in the case of a wobbling grain, hence we get 

^peak(T OT oo) « ^ZVak(T OT ~+ 0) « 1-34 Z/ peak (T CT 0). (80) 

The total power radiated scales as the fourth power of the angular momentum and at equal 
angular momentum it is ~ ten as large in the case of a wobbling grain, and therefore the 
total power radiated is roughly twice as large in the case of wobbling grains. 

This heuristic argument is in excellent agreement with results from detailed calculations. 
We show in Fig. [3] the difference in emissivity in the WIM environment. The peak frequency 
is enhanced by a factor of 1.33 and the total power by a factor of 1.9 for wobbling grains. 
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Figure 3: Effect of wobbling of axisymmetric grains on the spinning dust emissivity in the 
WIM environment. The spectra were produced with SpDust. 



Hoang et al [TU] studied the intermediate case where the internal relaxation temperature T ro 
is set to a finite value, and obtain similar results as they vary it from low values to large 
values. A similar heuristic argument could be made to estimate the peak frequency and total 
radiated power as a function of T^. 

Hoang et al [TTJ also studied the effect of triaxiality. They found an additional enhance- 
ment of the peak frequency and total power by up to the same factors (~ 30% and 2, 
respectively) for a large internal relaxation temperature and highly elliptical grains. 

One can therefore not neglect the fact that small PAHs are likely to be somewhat triaxial. 
The difficulty in properly accounting for this is that the exact distribution of ellipticities is 
largely unknown. 



6 Concluding remarks 

In this article we have reviewed the current status of spinning dust modeling, and tried to 
summarize the recent advances in this field since the seminal papers of Draine and Lazar- 
ian [21 [221 QUI [2H ITT] . In addition to refined calculations, the most important new effect 
accounted for recently is grain wobbling following frequent absorption of UV photons. The 
rotational dynamics of small grains of various shapes is now believed to be well understood, 
even if there remain uncertainties and simplifications in the implemented models. 

The accuracy of theoretical predictions remains mostly limited by our poor knowledge 
of the properties of small grains, namely their dipole moments, shapes and sizes, as well 
as their overall abundance, about which other observations give little information. This 
uncertainty can be turned into an asset, as one could potentially use the observed spinning 
dust emission (assuming it is the dominant AME process at tens of GHz frequencies) to 
constrain properties of small grains. 

Such a procedure can, however, only be accomplished if environmental parameters are 
very well known. Indeed, the gas density, temperature and ionization state as well as the 
ambient radiation field all affect the rotational distribution function of small grains in non- 
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trivial ways. In addition, the actual observable, the emissivity, depends upon the properties 
of the medium along the line of sight, and an accurate modeling of the spatial properties of the 
environment is also required. Unless the properties of the environment are well understood, 
it seems very difficult to extract dust grain parameters from observed spectra, due to the 
important degeneracies that are bound to be present for such a large parameter space. 

The view of the author is that significant advances in the field would be possible if 
several regions of the ISM were put under the scrutiny, not only of radio telescopes, but also 
of instruments at other wavelengths, in order to determine their detailed properties as much 
as possible and get rid of the uncertainties related to environmental dependencies. 

Finally, let us mention another potentially interesting avenue to probe the properties of 
emitting grains, namely the high-resolution spectral properties of the spinning dust spectrum. 
Indeed, even if the PAHs are classical rotators with large rotational quantum numbers, the 
line spacing remains relatively large for the smallest molecules (for coronene for example, ro- 
tational lines are spaced by about 0.33 GHz). A large number of different grains are probably 
present in the ISM, which results in a dense, quasi-smooth forest of lines. However, grains 
with a few tens of atoms might only be present in a limited number of stable configurations, 
or there might only be a fraction of possible grain configurations that lead to a significant 
electric dipole moment. If this were the case, radio observations with a narrow bandwidth 
should allow to detect some amount of bumpiness on top of a smooth spectrum. Even upper 
limits on the variability of the spectrum in the frequency domain should allow one to get 
some handle on the properties of small grains. A quantitative analysis of this issue will be 
the subject of future work. 
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